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Q ' The properties and the isospin dependence of the Hquid-gas phase transition in 



hot asymmetric nuclear matter have been investigated within the framework of the 
finite temperature Brueckner-Hartree-Fock approach extended to include the con- 
tribution of a microscopic three-body force. A typical Van der Waals structure has 
been observed in the calculated isotherms (of pressure) for symmetric nuclear matter 
implying the presence of the liquid-gas phase transition. The critical temperature of 
the phase transition is calculated and its dependence on the proton-to-neutron ratio 
is discussed. It is shown that the three-body force gives a repulsive contribution to 
^ . the nuclear equation of state and reduces appreciably the critical temperature and 

the mechanical instable region. At fixed temperature and density the pressure of 
\ asymmetric nuclear matter increases monotonically as a function of isospin asym- 

■ metry. In addition, it turns out that the domain of mechanical instability for hot 

asymmetric nuclear matter gradually shrinks with increasing asymmetry and tem- 
perature. We have compared our results with the predictions of other theoretical 
models especially the Dirac Brueckner approach. A possible explanation for the 
discrepancy between the values of the critical temperature predicted by the present 
non-relativistic Brueckner calculations including the three-body force and the rela- 
tivistic Dirac-Brueckner method is given. 
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I. INTRODUCTION 



The determination of the properties of nuclear matter as a function of tempera- 
ture, density and isospin asymmetry is one of the fundamental subjects in nuclear 
physics. For a long time it has been expected that, due to the Van der Waals nature 
of the nucleon-nucleon (NN) interaction, nuclear matter is likely to exhibit a liquid- 
gas phase transition at low densities and moderate temperatures^]. Phase transition 
is a phenomenon of great and extensive interests observed or expected in various 
fields, such as melting or boihng in condensed matters and atomic clusters^, 1^, 
decay of hot nucleiQQ and transition to qnark-glnon plasma in high energy heavy 
ion collisions and neutron stars p]. In heavy ion collisions at intermediate and high 
energies, hot asymmetric nuclear matter can be produced and during the expansion 
stage of the collisions the liquid-gas phase transition may occur depending on the 
temperatures and densities involved 3, 0]. Besides its general interest in nuclear 
physics, the equation of state (EOS) of hot asymmetric nuclear matter plays also 
an essential role in understanding many astrophysical phenomena related to the dy- 
namics of supernova explosions and the evolution of the protoneutron star formed 
in the latest stage of a type-II supernova collapse^l. It is therefore highly desir- 
able to predict on a microscopic basis the thermal properties of asymmetric nuclear 
matter at finite temperature which provide a starting point for the physics of heavy 
ions 3, [ill and protoneutron stars 12 1. 



Within various theoretical models such as the phenomenological non-relativistic 
Skyrme force models^, Q| and the relativistic mean field theory (RMT)^], the 
EOS of hot nuclear matter has been studied extensively. All of these investigations 
have predicted a Van der Waals behavior for infinite symmetric nuclear matter. The 
critical temperature Tc of the liquid-gas phase transition turns out to be in the range 
from 14 MeV to 20 MeV depending on the adopted models and NN interactions. 
Based on realistic NN interactions, the EOS of cold nuclear matter has been explored 
by many ph.*is. Within the v^iationalmethodQ , the non-relativistic Brueckner- 



Hartree-Fock (BHF) approach|l6l.ll7l| and the Dirac Brueckner (DB) theorvjl8llig|. 
However, the investigations of hot asymmetric nuclear matter from the microscopic 
models are relatively rare. In R.ef . |2fl| . the liquid-gas phase transition has been 
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predicted for symmetric nuclear matter by using the variational method based on 
the Argonne V14 two-body NN interaction plus a phenomenological three-body force 
and the critical temperature obtained is about 17.5 MeV close to the values from 
the Skyrme force calculations. On the contrary, it is reported in Ref. ^ from 
the finite temperature DB calculations an almost complete destruction of Van der 



•3, 



this unusual result in the DB 



Waals like behavior. As been discussed in Ref. 
framework has not been well understood and deserves deeper investigations. In 
a more recent workj23|, the Van der Waals structure has been obtained in the 
DB framework, but the predicted critical temperature is as low as ~ lOMeV. In 
Ref.Q, 

the properties of hot asymmetric nuclear matter has been studied in the 
finite temperature BHF approach by adopting the separable version of the Paris 
two-body force. As well known, within a non-relativistic microscopic framework 
three-body forces are decisive for reproducing the empirical saturation properties 
of cold nuclear matter and in the zero-temperature case their contribution shifts 
the EOS close to the one of the relativistic DB calculations In the present 
paper, our aim is to explore the EOS of hot asymmetric nuclear matter in the finite 
temperature BHF approach extended to include the contribution of the microscopic 
three-body force (TBF) constructed from the meson exchange current model 26[. 
We shall stress the isospin dependence of the liquid-gas phase transition and the 
effect of the TBF in particular its possible connection to the relativistic effect in 
the DB approach. In the present calculations, the Argonne Vis (AVis) potential {2?! 
has been adopted as the realistic two-body interaction. The paper is arranged as 
follows. In the next section we give a brief description of the theoretical approaches. 
The numerical results are presented and discussed in Section III. Finally a summary 
is given in Sect. IV. 

II. THEORETICAL MODELS 

A. Temperature-dependent EfTective Three-body Force 

The microscopic TBF adopted in t he p resent calculations is constructed from 
the meson-exchange current approach [2a| and its components are displayed dia- 



grammatically in Figd taken from Ref. 



. Four important mesons n, p, a and u 



are considered j28|. The TBF model contains the contribution of the two-meson ex- 
change part of the NN interaction medium-modified by the intermediate virtual ex- 
citation of nucleon resonances, the term associated to the non-linear meson-nucleon 
coupling required by the chiral symmetry, the simplest contribution rising from 
meson-meson interaction and finally, the two-meson exchange diagram with the vir- 
tual excitations of nucleon-antinucleon pairs. The meson masses in the TBF have 
been fixed at their physical values except for the virtual cr-meson mass which has 
been fixed at 540MeV according to Ref. This value has been checked to sat- 
isfactorily reproduce the AVis interaction from the one-boson-exchange potential 
(OBEP) modelQ- 

The other parameters of the TBF, i.e., the coupling constants 
and the form factors, have been determined from the OBEP model to meet the 
self-consistent requirement with the adopted AVig two-body force. The parameters 
of the TBF is given in Ref.j3]. A more detailed description of the model and the 
approximations can be found in Refs. 0, 2fil |. 

In the zero-temperature case, the TBF contribution has been included in the 
BHF calculations by constructing an effective two-body interaction via a suitable 



By extending 



average with respect to the third-nucleon degrees of freedom 
the standard scheme adopted for the zero-temperature case, one can reduce the TBF 
to a temperature dependent effective two-body force V^^{T) which reads in r-space, 

X (1 - ?7(r23,T))1^3(f/f2'f3'|fif2f3) (1) 

X 0„(r3)(l-r7(ri3,T))(l-77(r23,T)), 

where the trace is taken with respect to spin and isospin of the third nucleon. 
The function ri{r,T) is the defect function 0, Q, 13 which is defined as ri{r,T) = 
(f){r) — tp{r, T), where ip{r, T) is the correlated wave function for the relative motion 
of two nucleons in nuclear medium and 0(r) is the corresponding unperturbed one. 
A detailed description and justification of the above scheme can be found in Ref. 2^. 

It is worth stressing that the TBF itself, i.e. W3, is the same as the one adopted 
in our previous calculations for the zero-temperature case j 3] and it is independent 
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of temperature. However, in the finite-temperature case, tlie effective two-body 
force V^^{T) constructed from the TBF depends on temperature in an imphcitly 
comphcated way due to the medium effects. As a consequence the contribution of 
the TBF is expected to be more pronounced at finite temperature. It is clear from 
Eq.(0) that the temperature dependence of V^^{T) stems from the Fermi distribution 
f{k, p, P, T) and the defect function rj{r, T) which is strongly temperature dependent. 
We will return to this point in the next subsection. 



B. Finite Temperature Brueckner-Hartree-Fock Approach 

In general, three independent parameters are required to specify a given thermo- 
dynamical state of hot asymmetric nuclear matter, i.e., the total nucleon number 
density p, the isospin asymmetry parameter [3 = {pn — Pp)/p and the temperature 
T. At zero temperature, the neutron and proton Fermi momenta are related to their 
respective number densities p„ = {l + P)p/2 and pp = (1 — /9)p/2 by kp = [Stt^p,-]^/^ 
with T = p or T = n. The Brueckner-Bethe-Goldstone (BBG) approach for cold 
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i3- 



The extension to finite tem- 



asymmetric nuclear matter is described in Ref. 

perature is given in Ref.Q]. In the following, we give a brief review for completeness. 
The starting point of the BBG approach is the Brueckner reaction G matrix which 
satisfies the following generalized Bethe-Goldstone (BG) equationjlfi|. 

Gry {p,(3,T,Uj)=V + V}^ — — Gr,r' p, /?, T, Uj) , 2 

fcTte uj-er{ki)-er'{k2) 

where v = V2 + V^^{T) is the NN interaction and ut is the starting energy. In 
the present calculations the AVis potential is adopted as the bare two-body force 
V2 and the TBF contribution is included via the effective interaction V^^{T) given 
by Eq.((H). Since the defect function is determined by the G matrix, the effective 
force V^^ should be evaluated self-consistently with the BG equation at each step 
of the Brueckner iteration. In Eq.Q, the single particle (s.p.) energy is defined as 
e^(fc) = er{k, p, P,T) = ?i^k'^/(2m) + Urik, p, P,T). In the present calculations, we 
adopt the continuous choice |3fl| for the s.p. potential Ur{k) = Ur{k, p, P,T) since it 
is an natural choice for T 7^ 0| and it provides a much faster convergence of 
the hole-line expansion in the zero-temperature limit than the gap choice j33|. The 
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s.p. potential is calculated from the real part of on-shell G matrix, 

k' 

where the subscript A denotes the anti-symmetrization of the matrix elements. The 
finite temperature Pauli operator is simply an extension of the zero-temperature 



Q^,^ik,,k2) = Qr,r'{k,,k2,p,(3,T) = [l-f^{k^,p,(3,T)][l-Mk2,p,P,T)]. (4) 

In Eq.Q, the Pauli operator is applied only for the intermediate states of momenta 
ki,k2. The Fermi distribution for T 7^ is expressed as, 

-1 



fr{k,p,f3,T) 



1 + exp 



e-rjk) - Pr 

T 



(5) 



where pr = Prik, p, T) is the chemical potential. For any given density and tempera- 
ture, we can calculate the chemical potential pr from the following implicit equation 
self-consistently by iteration, 

(6) 



Pr = ^Efrik,P,(3,T) = W 

k k 



(er{k)-pr^ 

1 + exp - 



T 



The BG equation can be expended in partial waves as in the zero-temperature 
casejs^. The resulting equations in different partial waves are coupled with each 
other due to the angular dependence of the Pauli operator and the s.p. energy. To 
remove the angular dependence, we approximate the exact Pauli operator and s.p. 
energy in the BG equation by their angle-averaged values. For instance, the angle 
average of the Pauli operator is expressed as, 

{Qrr' (g, P,p,P,T)) = l r Sin 9d9 [1 - fr{k,, p, P, T)][l-frik2,p,P,T)], (7) 

2 Jo 

where q = {ki — k2)/2 and P = ki + k2 are the relative momentum and the total 
momentum of the two nucleons, respectively. 6 is the angle between q and P. In the 
case of T = 0, the angle-average of the Pauli o pera tor can be derived analytically 
due to the sharp Fermi distribution at T = ollfii IitI Whereas for T 7^ 0, 

the integral in Eq.© can not be worked out explicitly and one has to evaluate the 
angle-average of the Pauli operator numerically. In the present paper, the integral in 
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Eq.© is calculated numerically by the Gauss- Legendre method. We have checked 
that in the zero-temperature limit our numerical results can reproduce the analytical 
values with very high accuracy. 

The energy per nucleon of asymmetric nuclear matter at the BHF level of ap- 
proximation is given by, 

E{p,(3,T) = Y^Y.frik,P,P,T)^ + WY.Mk,p,(3,T)fAk',p,f3,T) 

T k t,t' k,k' 

X {kk' I Gr,r'{p, /9, T, e{k) + e{k')) \ kk')A. (8) 
At the mean field approximation, the total entropy S is expressed as|24| 
S = -J2 Eifrik, P, /3, T) In Mk, p, (3, T) + [1 - /.(fc, p, /3, T)] ln[l - Mk, p, /3, T)]} 

T k 

(9) 

The free energy F is calculated according to the standard thermodynamic relation 
F = E — TS and the pressure P is then extracted from the free energy by performing 
a numerical derivative, i.e., 

„ '■'■(l)„ 

According to Ref. |8| the criteria for mechanical stability can be straightforwardly 
expressed as follows 

The set of coupled equations (j2 1, (El) and ^ is referred to as the finite temperature 
BHF approximation (FTBHF) 24|. When the TBF contribution is included, this 
set of equations has to be solved self-consistently along with Eq.Q to get the G 
matrix. Hereafter we will call the FTBHF including the TBF contribution as the 
extended FTBHF. 

As a check of the extended model, we have calculated the speed of sound in cold 

/,p\l/2 

symmetric nuclear matter, i.e., ■^/c = ( , where c represents the speed of light 

and e refers to the total energy density. The results are shown in Tab. 1 for several 
densities. In the table the result including the TBF contribution is denoted by 
BHF(AVi8+ TBF), while the one without the TBF by BYiY{AVi^). It is seen that 
in both cases with and without the TBF, the calculated speed of sound fulfills the 
causality condition s/c < 1 in the relevant density region. Below a certain value of 
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density (~ 0.14fm ^ in the case without the TBF and ~ 0.11 fm ^ with the TBF) 
the speed of sound sjc becomes imaginary due to the mechanical instabihty. 

TABLE I: Speed of sound s/c in symmetric nuclear matter for several values of density in 
both cases, with and without the TBF. 



p(fm ^) 


0.15 


0.20 


0.25 


0.30 


0.35 


0.40 


0.45 


BHF(^yi8) 


0.01 


0.08 


0.15 


0.20 


0.25 


0.29 


0.35 


BHF(Ayi8+ TBF) 


0.07 


0.16 


0.24 


0.33 


0.42 


0.51 


0.59 



III. NUMERICAL RESULTS AND DISCUSSIONS 

The EOS of asymmetric nuclear matter is reported in Fig|2 for four values of 
asymmetry parameter /5 = 0,0.3,0.5,0.8, where the solid and dashed isothermal 
curves of pressure ( corresponding to T = 0, 8, 10, 12, 14, 16MeV from the bottom 
to the top) indicate the results with and without the TBF contribution, respectively. 
From the figure one can see that for both cases with and without the TBF, the EOS 
of symmetric nuclear matter (/3 = 0) displays a typical Van der Waals behavior, 
implying that the infinite nuclear system may undergo a liquid-gas phase transition. 
The critical temperature and density of the phase transition is determined by the 
condition (f^)^^ = (^;?)^^ = 0. From the calculated isotherms of pressure, one 
can extract the critical temperature Tc and the critical density pc- In the case with- 
out the TBF contribution, the obtained Tc is approximately 16 MeV which is in the 
range of 15 — 20 MeV predicted by the Skyrme-Hartree-Fock (SHE) calculations^^. 
As expected, the TBF gives a repulsive contribution to the nuclear EOS. In the 
zero-temperature case, the additional repulsion from the TBF improves greatly the 
predicted saturation density of cold symmetric nuclear matter. The equilibrium 
properties of cold nuclear matter have been reported in Ref . j^l • At higher temper- 
ature, the TBF effect becomes more pronounced due to the temperature dependence 
of its contribution (i.e., the effective force V^^{T)). Inclusion of the TBF contribu- 
tion reduces the critical temperature from ~ 16 MeV to ~ 13 MeV which is close 
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to the value of ~ 14 MeV from the relativistic mean field theory|8[. The reduction 
of Tc due to the TBF may be attributed largely to the medium and temperature 
dependence of the reduced effective interaction V^^{T), i.e., at a fixed temperature 
the TBF contribution is stronger for larger density and for a given density it is more 
pronounced at higher temperature. The calculated critical density pc of the phase 
transition is roughly 0.065 fm~'^ and 0.08 fm~'^ in the two cases with and without 



the TBF, respectively. Both va 
ment with other investigations 



ues are in the density range of po/3 ~ Po/2 in agree- 



13 



14 



j. The TBF reduces the critical density 
since its contribution to the nuclear EOS is repulsive and increases as a function of 
density. 

Now let us compare our result with the predictions of other approaches. In 
Ref.Q, 

the liquid-gas phase transition of hot nuclear matter has been studied 
with the variational method by adopting the Argonne V14 two-body NN interaction 
supplemented with a phenomenological TBF. The obtained value of Tc in Ref. [2^ is 
about 17.5 MeV which is higher than the present value Tc ~ 13MeV from the BHF 
approach when the TBF contribution is included. This disagreement concerning 
the critical temperature might be attributed to the two different approaches and 
to the different NN interactions adopted. As it has been noticed in Ref. 

Q, the 

density dependence of the symmetry energy predicted from the two methods differs 
remarkably from each other even when the same NN interaction has been used. The 
reason for the discrepancy between the results of the two methods remains largely 



unclear as discussed in Ref. 



In the relativistic DB framework, the EOS of hot nuclear matter has been inves- 
tigated in Refs. 



21 



23| where the extension to finite temperature has been achieved 



by the inclusion of the finite temperature Green's functions 



3J]. Such an extension 



coincides in the non-relativistic limit with the one adopted in the present FTBHF 



calculations. In Ref. 



2l|, the reported EOS does not display any Van der Waals 



like behavior in contrast to the results obtained from other approaches such as the 
non-relativistic BHF, the SHE, the variational method, and the RMT. As discussed 
in Ref. [2^ this discrepancy is not easy to be understood since the inclusion of the 
TBF in the BHF calculations is expected to give a similar behavior of the EOS to 
that of the DB approach. One possible reason for that is that in the calculations of 



10 



Ref. 



21[ , the A-resonance has been included exphcitly which produces an additional 



the 



T, ~ 10 MeV in Ref.|22| is considerably 



repulsion to the EOS and makes nuclear matter less bound. In a later work 
problem has been re-visited in the DB approach and a clear Van der Waals struc- 
ture is observed in the calculated EOS. However, the obtained critical temperature 

ower than the value of ~ 13 MeV in the 
22l 3511 that the main relativistic effect in the 



present work. As pointed out in Refs. 
DB approach is associated to the contribution (of the TBF) due to the 2(T-exchange 
process coupled to the virtual excitation of a nucleon-antinucleon pair which is re- 
ferred to as 2a-NN TBF in Refs.j^, l2^- In our previous investigations [3], it is 
shown that the most important relativistic correction to the nuclear EOS can be 
fairly well reproduced by the BHF calculation including only the 2a-NN contribu- 
tion of the TBF. This implies that the full TBF in the present calculations is not 
completely equivalent to the relativistic effect in DB approach especially in the finite- 
temperature case for which both the relativistic and the TBF effects are expected to 
become more pronounced as compared to the zero-temperature case. In order to get 
a deeper insight into the disagreement of the predicted critical temperature between 
the DB approach and the FTBHF approach with the TBF contribution, we separate 
the 2a-NN contribution from the full TBF. The results are given in FigOlwhere the 
solid curves are obtained by including only the 2a-NN contribution and the dashed 
curves by adopting the full TBF. It is seen that the contribution of the full TBF is 
somewhat less repulsive than that of the 2a-NN component. As a consequence, the 
value of Tc is lowered from ~ 13 MeV to ~ 11 MeV if only the 2a-NN contribution 
in the TBF is taken into account. This result is compatible with our previous one 
for the zero-temperature case 2^ and provides a possible explanation for the above 
mentioned discrepancy. 

As well known, pure neutron matter is unbound, therefore it is interesting to 
discuss the isospin dependence of the phase transition. By comparing the results 
for different isospin asymmetries (Fig|2K,b,c and d), one can see that as the isospin 
asymmetry increases, the Van der Waals structure of the EOS becomes less pro- 
nounced and the mechanical instable region where the pressure decreases as a func- 
tion of density, gets smaller. This indicates that the critical temperature drops down 
monotonically as the asymmetry increases in both cases with and without the TBF, 
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in agreement with the results obtained from the isospin lattice gas model js^]- The 
TBF gives a repulsive contribution to the EOS. The repulsion of the TBF becomes 
stronger as the density increases and consequently reduces the critical temperature 
in the whole range of asymmetry. When the asymmetry is high enough, the region of 
mechanical instability vanishes for all values of temperature. The extracted critical 
asymmetry (3^ for the disappearance of the mechanical instability is about 0.85 in 
the case without the TBF contribution. Inclusion of the TBF reduses (3c to about 
0.75. It is also seen that at a fixed asymmetry, the mechanical instable domain of 
the system gets smaller as the temperature becomes higher. 

In FigH is shown the pressure as a function of density at a fixed temperature 
T = 5 MeV for different values of asymmetry (3 = 0.0,0.2,0.4,0.6,0.8,1.0. In the 
figure, the solid curves and dashed curves correspond to the results with and without 
the TBF contribution. It can be observed from the figure that the pressure increases 
monotonically as the isospin asymmetry increases in agreement with the results 
of the RMT calculations p| where it is shown that for asymmetric nuclear matter 
{(3 7^ 0) which is a two-component system, the pressure cannot remain a constant 
during the liquid-gas phase transition. In the cases of (3 = 0.0, 0.2, 0.4, 0.6, the 
system exhibit a mechanical instability where the pressure is a decreasing function 
of density. However for high enough asymmetry, especially for pure neutron matter 
(/5 = 1), the pressure becomes an monotonically increasing function of density, which 
indicates the disappearing of the mechanical instability at all densities. The TBF 
gives a repulsive contribution in the whole asymmetry range < (3 < 1. 

In order to illustrate more clearly the TBF effect on the properties of the phase 
transition in asymmetric nuclear matter, we plot in FiglHlthe domain of mechanical 
instability in density-asymmetry plane for the fixed temperature T = 5 MeV. In the 
figure, the solid and dashed curves are obtained by adopting the AVis plus the TBF 
and the pure AVis two-body force, respectively. From the figure we can see that 
inclusion of the TBF contribution suppresses considerably the region of mechanical 
instability in isothermal asymmetric nuclear matter. In addition, the TBF effect 
on the lower-density boundary of the mechanical instable region is very small but 
it leads to a remarkable reduction of the upper-density boundary. This is expected 
since the TBF repulsion becomes stronger at higher densities. 
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IV. SUMMARY 



In the present work, we have introduced the microscopic TBF into the framework 
of the finite temperature BHF approach. Employing the extended theoretical model, 
we have investigated the EOS and the properties of the liquid-gas phase transition for 
hot asymmetric nuclear matter. For symmetric nuclear matter, the calculated EOS 
exhibits a clear Van der Waals behavior in the pressure-density plane, which implies 
the presence of a liquid-gas phase transition. The extracted critical temperature 
is about 16 MeV in the case without the TBF contribution. It is shown that the 
TBF gives a repulsive contribution to the EOS and its effect becomes stronger as 
increasing density. When the TBF contribution is included, the critical temperature 
turns out to be reduced by about 3 MeV from ~ 16 MeV to Tc ~ 13 MeV as 
compared to the corresponding two-body force predictions. The predicted value of 
Tc ^ 13 MeV is close to the value ~ 14 MeV obtained from the RMT^J] but it is 
appreciably larger than that of the relativistic DB approach. If including only the 
2cr — NN contribution of the TBF, we obtain a lower value of Tc ~ 11 MeV which is 
close to the value of about 10 MeV predicted by the DB calculations in Ref. [2^. This 
result is desirable since it may provide a possible explanation for the discrepancy 
between the critical temperature values obtained from the DB approach and the 
present BHF plus TBF calculations, i.e, except the 2ct — NN contribution, the other 
components of the TBF do not completely cancel among each others at least for 
the case of finite temperature and their net effect may becomes more pronounced at 
higher temperature. The critical density pc is found to be roughly 0.065 fm~^ and 
0.08 fm~^ in the cases with and without the TBF contribution, respectively. Both 
values are in the range between po/3 and po/2. 

For asymmetric nuclear matter, the isospin dependence of the properties of the 
liquid-gas phase transition has been studied. It turns out that the critical temper- 
ature decreases and the mechanical instable region gradually shrinks as the isospin 
asymmetry (3 increases in both cases with and without the TBF. The contribution 
of the TBF to the EOS is repulsive in the whole range of asymmetry < /5 < 1 
and results in a reduction of the critical temperature for the phase transition in hot 
asymmetric nuclear matter. Above a critical value of asymmetry, the mechanical 



13 



instability disappears for all values of temperature. The addition repulsion due to 
the TBF lowers the critical asymmetry for the disappearance of the mechanical in- 
stability and makes the asymmetric nuclear matter easier to be gasified. At fixed 

temperature and density the pressure in hot asymmetric nuclear system turns out 
to be an monotonically increasing function of asymmetry parameter. In addition, 
it is also shown that at a fixed asymmetry the mechanical instable region of hot 
asymmetric nuclear matter shrinks as increasing temperature. 
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FIG. 1: Diagrams of the microscopic TBF adopted for the present calculations, taken from 
Ref.[26] 

FIG. 2: Pressure as a function of density with six isotherms corresponding to T = 
0, 8, 10, 12, 14, IQMeV from the bottom to the top for different isospin asymmetries 
P = 0, 0.3, 0.5 and 0.8. The solid and dashed curves represent the results with and without 
TBF contribution. 

FIG. 3: Pressure as a function of density for symmetric nuclear matter at six values of 
temperature T = 0, 8, 10, 12, 14 MeV from the bottom to the top. The dashed curves are 
obtained by adopting the full TBF, while the solid ones by including only the 2a — NN 
contribution of the TBF. 

FIG. 4: Pressure as a function of density at a fixed temperature T = 5MeV for different 
asymmetric parameters P = 0,0.2,0.4,0.6,0.81. from the bottom to the top. The solid 
and dashed curves are the results with and without the TBF contribution, respectively. 

FIG. 5: Region of mechanical instability at a fixed temperature T = bMeV for both cases 
including the TBF (solid curve) and without the TBF (dashed curve) . 
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